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Abstract 

o 

■ We establish the existence of stationary Gibbsian point processes for interactions that 

act on hyperedges between the points. For example, such interactions can depend on 
Delaunay edges or triangles, cliques of Voronoi cells or clusters of fc-nearest neighbors. 
The classical case of pair interactions is also included. The basic tools are an entropy 
bound and stationarity. 
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1 Introduction 



Recent developments in statistical physics, stochastic geometry and spatial statistics in- 
volve Gibbs point processes with interactions depending on the local geometry of config- 
urations in W d . A prominent class of such interactions is based on the nearest-neighbor 
graph coming from the Delaunay triangulation. In biology, such systems are used to 
model interacting cells in tissues or foams |llj . |20j . In spatial statistics and stochastic 
. geometry, structured point patterns, point processes lying along fibers or regular Delau- 

nay tessellations have been studied via geometric Gibbs modifications; see pQ, [6] and |21| . 
A probabilistic motivation comes from stationary renewal processes: Since these can be 
characterized as Gibbs processes for interactions between nearest-neighbor pairs of points 
\16\ Section 6], Gibbs processes on M. d with Delaunay tile interaction can be viewed as a 
multi-dimensional counterpart [5]. 
^ ■ In this paper we consider general geometry-dependent interactions that are defined on 

a hypergraph structure £. For every point configuration oj in M. d , £(uj) denotes a set of 
hyperedges on uj, and the formal Hamiltonian of oo is given by 
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with a potential (p(rj, oj) which only depends on rj and the points of oj in some neighborhood 
of r\. This locality of <p will be called the finite horizon property, see (|2.ip below. This 
general setting includes all the above-mentioned cases and also the classical many-body 
interactions of finite range that are familiar from statistical mechanics [22 j. 

There is a principal difference between the geometric interactions considered here and 
the classical many-body interactions. Namely, suppose a particle configuration oj is aug- 
mented by a new particle at x. In the case of a many-body interaction, this is only 
influenced by an additional interaction term between x and the particles of lj, and the 
interaction between the particles of oj is not affected by x. In other words, the classical 
many-body interactions are additive. In our setting, the new particle at x typically al- 
ters the hyperedges around x completely: some hyperedges are created and some others 
are destroyed. This means that both £{oj U {x}) \ £(oj) and £{oj) \ £{oj U {x}) are non- 
empty, so that H(uj U {x}) and H(uS) each contain terms that are not present in the other. 
This phenomenon blurs the usual distinction between attractive and repulsive interactions. 
Moreover, if the potential ip(r],oj) is allowed to take the value oo (hard-exclusion case), we 
arrive at the so-called non-hereditary situation that a configuration oj is excluded although 
oj U {x} is possible. This last case makes it difficult to use an infinitesimal characteriza- 
tion of Gibbs measures in terms of their Campbell measures and Papangelou intensities. 
Nevertheless, such an infinitesimal approach was first used to prove the existence of Gibbs 
measures for Delaunay interactions by requiring geometric constraints on the interaction 
[H EJ HI |6] . quite different global approach, first used in [2] and based on stationarity 
and thermodynamic quantities such as pressure and free energy density, recently allowed 
to prove the existence of Gibbsian Delaunay tessellations for general bounded interactions 
without any geometric restrictions [8], and a similar approach could also be applied to 
quermass- interaction processes 0. 

In this paper we address the existence problem for the general formalism of hypergraph 
interactions introduced here. Our approach is global as in [3 E] and leads to a significant 
improvement of the existing results. In particular, we establish the existence of Gibbsian 
Delaunay tessellations for non-bounded and hard-exclusion potentials. In the classical 
context of stable many-body interaction of finite range, our results permit to relax the 
superstability assumption. Basic ingredients of the proof are an entropy bound to exploit 
the compactness of the level sets of the entropy density, and a somewhat delicate control 
of the range of the interaction, which takes advantage of stationarity. 

The general setting of Gibbs measures for hypergraph interactions is introduced in 
Section 2. Section 3 contains our assumptions and the existence theorems. Section 4 
offers a series of examples that includes many-body interactions of finite range as well as 
interactions on Delaunay tiles and Voronoi cells and between /c-nearest neighbors. The 
proofs of the main results follow in Section 5, and an appendix is devoted to measurability 
questions. 

2 Preliminaries 

2.1 Point configurations and hyperedge interactions 

Consider the Euclidean space R rf of arbitrary dimension d > 1. Subregions of M. d will 
typically be denoted by A or A and will always be assumed to be Borel with positive 
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Lebesgue measure |A| resp. |A|. We write A (e R if A is bounded. A configuration is a 
subset oj of R d which is locally finite, in that wnA has finite cardinality Na{oj) = #(wnA) 
for all Aglf The space fl of all configurations is equipped with the cr-algebra J- that 
is generated by the counting variables N A with A <e R d . It will often be convenient 
to write wa in place of oj D A. As usual, we take as reference measure on (J), J 7 ) the 
Poisson point process IF of an arbitrary intensity z > 0. Recall that IF is the unique 
probability measure on (fl,T) such that the following holds for all A <s R d : (i) N A is 
Poisson distributed with parameter z\A\, and (ii) conditional on N A = n, the n points in 
A are independent with uniform distribution on A, for each integer n > 1. 

Next, let fl f C fl denote the set of all finite configurations oj (which means that 
< oo), and J 7 / = F\n f the trace cr-algebra of T on flf. The product space $7/ X fl 
carries the product cr-algebra T$®T. For each A C R d we write fl\ = {oj G : oj C A} for 
the set of all configurations in A, pr^ : oj — > oj\ = u D A for the projection from fl to fl\, 

= F\n A for the trace cr-algebra of T on fl\, and F\ = pr^ 1 ^ C T for the cr-algebra 
of all events that happen in A only. The reference measure on (fl\, F'^) is 11^ := IF opr^ 1 . 
Finally, let = ('&x) x eR d be the shift group, where "& x : fl — >■ fl is the translation by the 
vector -x G R d . By definition, N A (tf x uj) = N a+x (lo) for all A m R d . 

The interaction of points to be considered in this paper will depend on the geometry 
of their location. This geometry will be described in terms of a hypergraph, and the 
interaction potential will be defined on the hyperedges. 

Definition. 

• A hypergraph structure is a measurable subset £ of fl f x fl such that rj C lo for all 
(rj,uj) G £. If (rj, oj) G £ , we say that 17 is a hyperedge of w, and we write r] G 

• A hyperedge potential is a measurable function if from a hypergraph structure £ to 
1U {oo}. 

• A hyperedge potential if (or, more explicitly, the pair {£ , 99)) is called shift-invariant 
if 

($ x r],'& x cu) G £ and <p('& x ri,'& x uj) = <p(r],uj) for all (r/,o;) G £ and x G M d . 

• Let us say that ip (or the pair (£, </?)) satisfies the /imte horizon property if for each 
(77, w) G £ there exists some A <s R d such that 

(2.1) (f]^) ^ ^ an d <p(rj,uj) = 99(77,0;) when a) = uj on A. 

W^e wi// assume throughout this paper that the hyperedge potential ip under consideration 
is shift- invariant and exhibits the finite horizon property. Moreover, for notational conve- 
nience we set if = on £ c . Since £ is measurable, this does not affect the measurability 

Of if. 

The domain £ of ip can be considered as a rule that turns each configuration oj into 
a hypergraph (oj,£(oj)). Both £ and ip are not affected by translations. Moreover, the 
presence of a hyperedge r/ G £(oj) and the value of <p(r], oj) can be determined by looking at 
w in a (sufficiently large but) bounded neighborhood A of rj, called the horizon of r/ in oj, 
which in general depends on both rj and oj. Note that in general there is no minimal such 
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horizon. To obtain a standard choice of A one can take the closed ball = B(g v , r VjU ) 
with center at the gravicenter g v of 77 and radius r^ )W chosen smallest possible. Finally, 
we note that the concept of hypergraph structure is similar to that of a cluster property 
as introduced in [25]. Here are two examples the reader might keep in mind. Further 
examples will follow in Section SJ 

Example 2.1. Many-body interactions of bounded range. Let r > and 

LC r = {(77,0;) : 77 C oj, diam(77) < r, oj £ fi} 

be the locally complete graph. Thus, for each oj £ Q, LC r (w) consists of all hyperedges 
between points of distance at most r. If we assume that (p(rj,un) only depends on 77, 
we are in the classical situation of many-body interactions of range r. The finite horizon 
property for (77,0;) £ LC r then holds for arbitrary A D 77. If LC r is restricted to hyperedges 
of cardinality two, we arrive at the familiar pair interactions of statistical mechanics. 

Example 2.2. Delaunay potentials. The set Del of Delaunay hyperedges consists of all pairs 
(ij,oj) with 77 C oj for which there exists an open ball B(rj,uj) with dB(r], oj) n u = 77 that 
contains no points of oj. For k = 1, . . . , d+1 we write Delfc = {(77, oj) £ Del : #77 = k} for 
the set of all Delaunay simplices with k vertices. Clearly, Del and Del^ are hypergraph 
structures. It is possible that the convex hull of a set 77 £ Del (a;) is not a simplex, namely 
when 7/ consists of four or more points on a sphere with no point inside. However, this 
is an exceptional case, which occurs only with probability zero for our Poisson reference 
measure Tl z . Note that B(tj,uj) is only uniquely determined when #77 = d+1 and 77 is 
affinely independent. 

The simplest class of Delaunay hyperedge potentials consists of pair interactions of the 
form (p(rj,oj) = (p(\x — y\) for r/ = {x,y} £ Del2(w). Such a <p satisfies the finite horizon 
property (|2.ip with A = B(t],lu) for any ball B(t],lu) as above. 

An example of a potential <p(ri, oj) on Del2 which does not only depend on 77 but also on 
oj is (p(j],oj) = 0(Vor w (x), Vor w (y)). Here we write Vor w (x) for the Voronoi cell associated 
to a point x £ oj, viz. the set 

(2.2) Vor w (a;) := {y £ R d : \x - y\ < \x - y\ View} 

of all points of M d which are closer to x than to all other points of oj. It is well-known 
that the Voronoi cells form a tessellation |24j . Also, any two points of a configuration 
are connected by a Delaunay edge if and only if their Voronoi cells have a non-trivial 
intersection. That is, 

(2.3) {x, y} £ De\ 2 (oj) #(Vor w (x) n Vor w (y)) > 1. 

This reveals that the Delaunay graph is a nearest-neighbor graph. The potential above 
satisfies the finite horizon property f)2. 1 1) with A equal to the closure of the set 
UggDei(w) £nr;^0 oj), provided the cells Vor^(x) and \/or UJ {y) are bounded. The pro- 
viso is necessary because unbounded Voronoi cells, which can occur at the "boundary" of 
oj, are not protected by the points in A. So we must exclude from De^O*;) all edges {x, y} 
for which Vor w (x) U Vor w (y) is unbounded. 
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2.2 Gibbs measures for hyperedge potentials 

Our objective here is to introduce the concept of a Gibbsian point process for an activity 
z > and a hyperedge potential ip defined on a hypergraph structure £. First we will 
introduce the Hamiltonian for a bounded region A (<= M. d with configurational boundary 
condition uj G £1. This requires to consider the set of hyperedges rj in a configuration uj 
for which either rj itself or p(rj, uj) depends on the points of to in A. Specifically, we set 

(2.4) £\(oj) = {rj G £(oj) : ^(r/, ( U wa c ) / fiv^) f° r some ( G ^a}. 

Recall the convention that ^ = 0on £ c . So, if r\ G £(w) \ £(C U wa c ) for some ( £ Oa then 
either r\ £ £a(u), or rj is irrelevant for w, in that ip(rj,uj) = 0. The Hamiltonian in A wii/i 
boundary condition uj is then given by the formula 

(2.5) H A>UJ (0-= Yl v(^CUwac) forCG^A, 

r?£fA(CUw A c) 

provided this sum is well-defined. As usual, we also consider the associated partition 
function 

To ensure that these quantities are well-defined, we impose the following condition on the 
boundary condition uj. Let uo~ = (—(f) V be the negative part of <p. 

Definition. A configuration uj G is called admissible for a region A (s M. d and an activity 
z > if 

H\ w(0 := XT ¥' _ ( r ?' C U W A C ) < °o for n A -almost all ( G SIa 

(so that -£/a,w is almost surely well-defined), and < Z z k w < oo. We write Q^' z for the set 
of all these uj. 

We note that, in contrast to the standard setting of statistical mechanics, the partition 
function is not automatically positive because, in the present setting, H A ,w{^) is not nec- 

A z 

essarily finite. For , we can define the Gibbs distribution for (£, ip, z) in a region 

A (s W 1 with boundary condition uj as usual by 

(2.6) GIJF) = [ 1 F (C U u A c) e-^(C) n A (d() / , 
where F G T is arbitrary. 

Definition. Let £ be a hypergraph structure, a hyperedge potential, and z > an 
activity. A probability measure P on _F) is called a Gi66s measure for £, <p and z if 
P(nt ,Z ) = 1 and 




for every A <s R d and every measurable / : fi — »■ [0, oo[. 
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The equations (|2.7p are known as the DLR equations (after Dobrushin, Lanford and 
Ruelle). They express that G^ uJ (F) is a version of the conditional probability P( F \F\c)(uj). 
We will be particularly interested in Gibbs measures that are stationary, that is, invariant 
under the shift group = ( r d x ) xG -^d. We write &q for the set of all 0-invariant probability 
measures P on (fl,F) with finite intensity i(P) = J N^^ddP, and @Q(ip,z) for the set 
of all Gibbs measures for ip and z that belong to &q. We conclude this section with a 
discussion of measur ability questions. 

Remark 2.1. Measurability. The quantities introduced above are not measurable with 
respect to the underlying cr-algebras defined so far, but only with respect to their universal 
completion. Specifically, for each cx-algebra A let A* be the associated cr-algebra of all 
universally measurable sets, i.e., of the sets which belong to the P-completion of A for all 
probability measures P on A; see [5j pp. 36 & 280]. It is then the case that, for each A <g M. d , 
the Hamiltonian (Ci^) l— #Aw(C) i s measurable with respect to (J 7 ^ (g) J~a c )*- Likewise, 
the partition function lj — > is measurable with respect to F^ c , and Vlt' z G J 7 ^. 

Moreover, (oo,F) — > G\ U1 {F) is a probability kernel from (O*' 2 , ^Xc| n A,«) to (SI, J 7 ). All 
this will be proved in the appendix. We will therefore identify all probability measures in 
this paper with their respective complete extension. This convention underlies already the 
preceding definition of a Gibbs measure. 

3 Hypotheses and results 

Having defined the concept of Gibbs measure for a hyperedge potential we now turn to our 
main theme, the existence of such Gibbs measures. Let us state the conditions we need. 
In the subsequent section we will provide a series of examples for which these conditions 
are met. 

We begin with an assumption which says that hyper edges with a large horizon require 
the existence of a large ball with only a few points. This will imply that the Hamiltonian 
H\,lo depends only on the points of oj in a bounded region <9A(w), and can be viewed as 
a sharpening of the finite horizon property (|2.ip . 

(R) The range condition. There exist constants £r, ur G N and 5r < oo such that for 
all (77, oj) € £ one can find a horizon A as in f)2. lj) satisfying the following: 
For every x, y G A, there exist £ open balls Bi,...,Bg (with t < £r) such that 

- the set L)| =1 jB( is connected and contains x and y, and 

- for each i, either diami?j < 5r or Nb^uj) < ur. 

Note that (R) is trivially satisfied when all horizon sets can be chosen to have uniformly 
bounded diameters. For instance, this holds in Example l2.11 The use of the range condition 
(R) will be revealed by Proposition 13.11 below, which states that the following finite range 
property holds almost surely for nondegenerate P G 

Definition. Let A d M. d be given. We say a configuration oj G ft confines the range of if 
from A if there exists a set dA(ui) <s M. d such that ip(i], £ U &a c ) = V 9 ( r ?) C U w A c ) whenever 
Cj = oj on dA(oj), £ G S7a and r} G £\{C U wa c )- In this case we write u G Sl^ r . 9A(w) is 
called the w-boundary of A, and we use the abbreviation d\u = ojgu^y 
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Given any uj G Q^., we assume in the following that dA(ui) = A r \ A, where A r is the 
closed r-neighborhood of A and r = r A ^ is chosen as small as possible. Moreover, for 
uj G f2„ we have 

(3.1) H A>UJ {()= <f{V,(Ud A uj), 

and this sum extends over a finite set. This means that the first assumption in the 
definition of admissibility for A is satisfied. Here is the proposition announced above. It 
will follow from Proposition 15.41 below. 

Proposition 3.1. Under (R), for each A <s= M d there exists a set G J- A c such that 
ft£ C fi£ and P(fi£) = 1 for all P G ^ e with P({0}) = 0. 

Our next assumption is stability, the standard assumption that ensures the finiteness of 
all partition functions. In our setting, a somewhat modified definition turns out to be 
suitable. 

(S) Stability. The hyperedge potential ip is called stable if there exists a constant cs > 
such that 

(3.2) H A>UJ (0>-cs#((Ud A uj) 

for all A <e R d , ( G U A and uj G 

In Remark 13 . 41 below we will show that this definition is a natural extension of the familiar 
concept of stability in statistical mechanics. Complementary to the lower bound provided 
by stability, we will also need a further condition that provides at least a partial upper 
bound for the Hamiltonians. This is because in the extreme case when ip is constantly 
equal to oo we have . = 0, so that the definition of Gibbs measures is meaningless. 
Let M G M. dxd be an invertible d x d matrix and consider for each k G 7L d the cell 

(3.3) C{k) := {Mx G R d : x - k G [-1/2, l/2[ d }. 

These cells together constitute a periodic partition of M. d into parallelotopes. For example, 
the columns Mi, ... , of M might form an orthogonal basis of M. d or, for d = 2, define 
the sides of an equilateral triangle. For brevity we write C = C(0). Let T be a measurable 
subset of Qc \ {0} and 

(3.4) r = G n : t? M fc(wc(fc)) G r for all k G Z d | 

the set of all configurations whose restriction to an arbitrary cell C(k), when shifted back 
to C, belongs to T. We call each uj G T pseudo-periodic. The required control of the 
Hamiltonian from above will then be achieved by the following assumption on the joint 
behavior of £ , (p and z. Recall that r AtU1 was defined before (|3.ip . 
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(U) Upper regularity. M and T can be chosen so that the following holds. 



(Ul) Uniform confinement: T C f2„ for all A <s M. d , and 



rr := sup sup < oo. 



(U2) Uniform summability: cp~ := sup ^~j77^a " < 00 > 

wer J7e£(w):?7ncV0 

where r) :={£;£ Z d : 77 n C(fe) 7^ 0} and is the positive part of 93. 

(U3) Strong non-rigidity: n^.(r) > e Cr , where cr is defined as in (U2) with 

ip in place of <p + . 

Hypothesis (Ul) states that the configurations in T confine the range of tp in a uniform 
way. So, for oj £ T, the w-boundary <9A(w) of a set A <g M d is contained in the rp-boundary 
<9 r A := A Tr \ A, and the cardinality of d^oj ist not larger than that of d\ui := ojqt^. On 
the other hand, condition (U2) provides a uniform upper bound for the local Hamiltonians 
.Ha,- on T. This bound is of order cr|A|, cf. (|5.8p below. Finally, hypothesis (U3) holds 
for all z in an interval of the form ]zo,oo[, provided that 11^ (r) > for some (and thus 
all) z > 0. Indeed, since ^ T it follows that 

(3.5) e ^ c i n£(r) = Y.T\ • • • / lr ^ Z1 ' • • • ' Xk ^ dxi --- dXk 

k=l ' ^ c •' c 

is then a strictly increasing function of z. We emphasize that condition (U) imposes 
conflicting demands on T. While (Ul) and (U2) suggest to choose the set T as small as 
possible, (U3) requires that T is not too small. The point will be to choose a set T that 
satisfies all requirements simultaneously. Here is our main existence theorem. 

Theorem 3.2. For every hypergraph structure E, hyperedge potential <p and activity z > 
satisfying (S) ; (R) and (U) there exists at least one Gibbs measure P £ &q((p,z). 

In some cases it is difficult to satisfy hypothesis (U3) when z is small. This occurs, 
for example, when a hard-exclusion hyperedge interaction enforces a minimal number of 
points per cell; see Proposition 14.21 But a slight variation of proof allows to establish the 
existence of a Gibbs measure for every z > also in this case. Assumptions (Ul) and 
(U3) are replaced by conditions (Ul) and (U3) as follows. 

(U) Alternative upper regularity. M and T can be chosen so that the following holds. 

(Ul) Lower density bound: There exist constants a, b > such that #(C) > a|A|— 6 
whenever £ £ £2/ is such that H\ tilJ (Q < 00 for some C C A <s R d and some 

uj £ r. 

(U2) = (U2) Uniform summability. 
(U3) Weak non-rigidity: Ii z c (T) > 0. 

Here is the modified existence theorem. 
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Theorem 3.3. A Gibbs measure P G &q(<p,z) exists also under the hypotheses (S), (R) 
and (U). 

It is often natural to choose T as the set of all configurations that consist of a single 
point in some Borel set A C C. So we define 

(3.6) T A = {( G n c : C = {x} for some x G A}. 

For r = T A , the assumptions (U) and (U) are respectively called (U A ) and (U A ) . In 
particular, (\J2 A ) and (\J3 A ) take the simpler form 

(U2 A )c+:=sup £ ^^<oo and 

(U3 A ) >e c ^. 
We then have the following corollary of Theorems 13.21 and 13.31 

Corollary 3.4. A Gibbs measure P G &@((p,z) exists under the hypotheses (S), (R) and 
either (U A ) or (U A ). 

We conclude this section with a series of comments on our assumptions and on the 
extension to marked particles. 

Remark 3.1. Bounded horizons. The conditions (R) and (U A ) hold as soon as c^({0}, {0}) 
is finite and {£ , <p) has bounded horizons, in that there exists some < oo such that 
i"rj,ui < r v for all (r],ui) G £. (The notation r V)U1 was introduced after (|2.ip .) Indeed, 
condition (R) holds trivially with 5r = 2r v . As for (U A ), let M = aE, where a > 2r^ 
and E is the identity matrix. Let A = B(0,b) be a centered ball of radius b < a/2 — r v . 
For r = T A , condition (Ul) holds with rr = r^. Moreover, by the choice of a and b, each 
7] G £{oj) with ui G T must be a singleton {x}, so that ip(r],ui) = ip({x},{x}). In view 
of the shift-invariance of <p, this means that (XJ2 A ) holds with cj[ = i/? + ({0},{0}) < oo. 
Finally, (U3" 4 ) holds if a and b are in fact chosen so large that also irzb 2 > e^ ^ ^ . 

Remark 3.2. Scale-invariant potentials. Suppose £ and cp are scale-invariant in the sense 
that (rn,rui) G £ and (p(rn,ru>) = (p(r),ui) for all (r],ui) G £ and r > 0. Here, rr] = {rx : 
x G rf\ and rw = {rx : x G ui}. Then Theorem 13.21 is still valid when assumption (U3) 
is replaced by (U3). Indeed, the scale invariance of (p implies that the image of a Gibbs 
measure for tp and z under the rescaling ui — > rui is a Gibbs measure for ip and zr~~ d . So, 
it is sufficient to have the existence of a Gibbs measure for large z, and this follows from 
the remark around (|3.5p . 

Remark 3.3. Stability via sublinearity of the hypergraph. We say that a hypergraph £ is 
sub linear if there exists a constant C < oo such that #<?(w) < C #(w) for every finite 
configuration ui. In this case, the stability is ensured by requiring that the hyperedge 
potential ip is bounded below, in that 

(3.7) <p(r],ui) > -c v 

for some c v < oo. If the sublinearity of the hypergraph structure fails, the stability can 
simply be achieved by requiring that the potential ip is nonnegative (i.e., = 0). For 
example, for d = 2 it follows from Euler's formula that the cardinalities of the Delaunay 
edges and triangles are sublinear [8] , so that the stability follows directly from (|3,7p . 
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Remark 3.4. Stability: comparison with the classical case. Consider the hypergraph struc- 
ture £ = LC r of Example EI] describing many-body interactions of range r. In contrast to 
Example 12 . 2 1 where the Delaunay tiles depend on the presence of further particles, it is then 
meaningful to define the energy of a finite configuration ( £ ttf by H(() = Ylr)€£(C) < ^^)" 
The classical stability condition asserts that H(() > —cs#(() for all C £ see [22J, for 
example. This follows from (S) by choosing A3( and u = 0. Conversely, the Hamilto- 
nian (|2.5|) is equivalent to the Hamiltonian Ha,w(C) := -^(C U d\uj) in the sense that the 
associated Gibbs distributions coincide (at least when ip < oo), and the classical stabil- 
ity assumption for H gives (S) for H^ tUJ . This shows that hypothesis (S) is essentially 
equivalent to the classical concept of stability. 

Remark 3.5. Sub-hypergraph potentials. Consider a shift-invariant sub-hypergraph struc- 
ture £' C £ of £ and a hyperedge potential p on £, and let <p' be its restriction to £'. In 
general, <p' does not satisfy the finite horizon property, but let us assume it does. Which 
of the assumptions (R), (S), (U) and (U) on p are inherited by p'? It is clear that 
assumptions (R), (Ul), (U2) and (U3) are hereditary. Assumption (U3) remains also 
valid for p', but for a different range of values of z because the constant cr is different 
in general. Assumption (S) is lost, but a positive exception is the case of Remark 13.31 
when stability follows from the sublinearity of £ and the lower boundedness of p; these 
properties are obviously inherited by p' . Assumption (Ul) is lost in general. 

Remark 3.6. Upper regularity in Delaunay models. For potentials acting on the Delaunay 
graph, the matrix M and the set T in hypotheses (V A ) and (U A ) will be chosen as follows. 
Let M be such that |Mj| = a > for i = 1, ... ,d and <(Mj, Mj) = n/3 for i ^ j, and 
let r = T A with A = B(0, b). If b < qqo for some sufficiently small constant qq > 0, the 
Delaunay neighborhood of a point x in a configuration w£T contains a minimal number 
of points denoted by 7^. For d = 2 one can take qq = \/3/6 and has 72 = 6, and the 
Delaunay neighborhood of the unique point xt in to n C{k) consists of the unique points 
xi in uinC(l) with I- he {(-1, 0), (-1, 1), (0, 1), (1, 0), (1, -1), (0, -1)}. For d > 2, 7d is 
less easy to determine but it is clearly not larger than 3 d — 1, the value corresponding to 
the case M = aE. 

Remark 3.7. Extension to the marked case. The preceding results can be easily extended 
to the case of particles with internal degrees of freedom, or marks. Let S be an arbitrary 
Polish space with Borel u-algebra S and reference probability measure fi. X serves as the 
space of marks. That is, each marked point is represented by a position x 6 M. d and a 
mark a € S, and each configuration w is a countable subset of M. d x S having a locally 
finite projection onto M. d . The role of the reference measure on the configuration space f2 
is taken over by the Poisson point process IF with intensity measure zX <8> /i, where A is 
Lebesgue measure on W d . The translations "Q x act only on the positions of the particles 
and leave their marks untouched. We do not discuss the further formal details here, which 
are standard and can be found in [19] or [15] . for example. What we want to emphasize 
here is that all definitions and results above carry over to this setting without any change, 
provided it is understood that all regions A or A in M d always refer to the positional part 
of a configuration. For example: the notation wa now stands for wn(AxS); a set A <s M. d 
is the horizon of (rj, u) £ £ if (77, cD) G £ and p(rj,u}) = p(i],uj) whenever Co = to on A x S; 
and the condition in Remark [3 . 1 1 should now read sup CTe 5 ^({(0, a), {(0, a)}) < 00 for some 
Borel set S C S with fj,(S) > 0. 
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4 Examples 



In this section we present a series of examples to illustrate our general existence results. 
These examples satisfy the assumptions of Theorem 13.21 or 13.31 and many of them have 
been introduced in practical or theoretical papers without justification or with justification 
in some partial cases. 

To sort the examples of this section we will distinguish whether or not the potential 
ip(r],u)) depends explicitly on u. If not, we speak of a pure hyperedge potential. Otherwise, 
the finite horizon property implies that (p(rj, u) actually depends only on some points of 
u) \ r) close to rj, which is expressed by speaking of a neighborhood- dependent hyperedge 
potential. But note that this distinction is merely a matter of how the interaction is 
represented. In the pure case, the extended potential (rj,uj) — > tg(r],uj) cp(rj) on £lf x 
clearly does depend on w, while in the other case one can often include the neighboring 
points of a hyperedge into an enlarged hyperedge to obtain a pure hyperedge potential. 

Most of the following examples are based on the Delaunay graph. For simplicity, we 
will then often confine ourselves to the case d = 2, in which the stability is ensured by 
Remark l3.3l as soon as ip is bounded from below. But the reader should note that analogous 
results hold also in higher dimensions when ip is nonnegative, so that (S) is trivial. 

4.1 Pure hyperedge interactions 

In this subsection we consider examples of hyperedge potentials ip which only depend on 
the first parameter, so that ip(rj,uj) = (p(rj, rj) =: <p(rj). 

4.1.1 Many-body interactions of finite range 

Let r > and £ = LC r . be the locally complete graph of Example 12. 1\ and suppose 
that ip(rj,io) = p{rf). Remark 13.11 then shows that a Gibbs measure exists as soon as the 
potential ip is stable and y({0}) < oo. By Remark 13. A\ the first condition is equivalent to 
the classical stability assumption, and the second is necessary for defining Gibbs measures. 
So, as was observed first in |144 Remark 4.2], the techniques used here allow to weaken the 
superstability assumption of Ruelle's classical existence result [23], provided the interaction 
has finite range. But our techniques neither allow to treat the case of infinite range nor 
to rederive Ruelle's probability estimates. 

An example of a stable but not superstable many-body interaction is the so-called 
quermass interaction. In space dimension d = 2, the associated Hamiltonian for a con- 
figuration w is a linear combination of area, perimeter and Euler-Poincare characteristic 
of the union of all discs of fixed radius that are centered at the points of uj. By the 
additivity of the Minkowski functionals, this Hamiltonian can be expressed in terms of a 
many-body interaction of finite range. Its stability is proved in [18] . but the superstability 
fails. The existence of Gibbs processes for this interaction has been proved in [7] by the 
same methods used here. 

4.1.2 Delaunay edge interactions 

Here we consider two classes of potentials cp on the hypergraph structure Del2 in M. 2 
which are bounded below and depend only on the length of the Delaunay edges, in that 
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(p({x,y}) := <f>(\x— y\) when {x, y} G Del2- Such potentials have been studied in [TJ [2j SIS]. 
We improve the existence results of these papers. 

4-1-2.1 Polynomially increasing edge interactions. Suppose that 

(4.1) <j)(£) < kq + Ki£ a for some constants kq > 0, k\ > and a > . 
Recall the definition of qq in Remark 13.61 

Proposition 4.1. Lei d = 2 and (p be a pure edge potential on De^ which is bounded 
below and satisfies (|4,ip . Then there exists at least one Gibbs measure for <p and every 
activity 

z > (l+2 eo )e 3KO (3ae 2 K 1 /2) 1 / a /(7TQ 2 ). 

Proof. We apply Corollary 13.41 with the assumption (\J A ). The stability condition (S) 
holds because of Remark 13.31 We know further from Example 12.21 that every edge r] G 
Del2(w) has the finite horizon B(r/,oj). This shows that assumption (R) is satisfied. 
Concerning assumption (\J A ), let M and T A be defined as in Remark 13.61 with b = g^a 
and a to be chosen later. The assumptions (U1 A ) and (\J2 A ) are then trivial. We also 
find that c A < 3(k + K ia a (l+2Q ) a ) . So, (U3 A ) holds as soon as z > K e Kiaa /a 2 , 
where Kq = e 3K ° / '(it Qq) and K\ = 3ki(1+2^o)°- Optimizing over a we find the value 
a = (2/ 'aK\) l / a , which gives the sufficient condition z > Ko(aKie 2 /2) l / a . □ □ 

The bound for z in Proposition 14.11 is certainly not the best possible. However, in the 
bounded case K\ = it reduces to the trivial condition z > 0. As for condition (|4.ip . at 
first sight it might seem surprising that <p(£) does not necessarily converge to zero when 
£ tends to infinity, but is even allowed to converge to oo. However, <p(£) should not be 
compared with a pairwise interaction that must decay over large distances, but with a 
nearest-neighbor interaction between oscillators that form an elastic lattice. A potential 
like the harmonic interaction cf)(£) = £ 2 is therefore quite natural. 

4-1.2.2 Long-edge exclusion. Suppose there are constants < £q < £\ < £2 such that 

(4.2) sup </>(£)< 00 and </>(£) = 00 if £ > £ 2 . 

Proposition 4.2. Let d = 2 and (p be a pure edge potential on De^ which is bounded 
below and satisfies (|4.2p . Then there exists at least one Gibbs measure for ip and every 
z>0. 

Proof. We apply Corollary 13 .41 with the assumption (V A ). Hypotheses (S) and (R) follow 
as in Proposition 14.11 Condition (\J1 A ) is satisfied since the long-edge exclusion condition 
ensures a minimal density of points. The values a and b in the definition of M and T are 
now chosen such that each {x, y} G Del2(w) (with uj G V ) satisfies £q < \x — y\ < £\. 
Then (U2 A ) follows, and (U3 A ) is obvious. □ □ 

Since 4>{£) may be equal to 00 for £ less than some r^ < £q, the present example includes 
the classical case of a hard-core interaction which forces the points to keep distance at 
least ro > from each other. 
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Remark 4.1. Let us consider whether or not Propositions l4.1l and l4.2l remain valid when the 
Delaunay graph is replaced by a subgraph, as was discussed in Remark 13.51 One possible 
subgraph is the Gabriel graph. By definition, the Gabriel edge set Gab2(w) consists of all 
edges {x,y} 6 Debfw) for which the open disc with center (x + y)/2 and radius \x — y\/2 
contains no point of ui. Since Gab2 is shift-invariant and local, the restriction of <p to Gab2 
remains a shift-invariant hyperedge potential satisfying the finite horizon property (|2.ip . 
By Remark 13.31 (S) is inherited when we pass from De^ to Gab2. In view of our choice 
of M and T, the constant cr remains unchanged, so that (U3) is also inherited without 
modification of the valid range for z. This means that Proposition 14.11 holds also for the 
Gabriel graph. Unfortunately, condition (Ul" 4 ) is lost, so that Proposition 14.21 does not 
carry over to the Gabriel graph. Another example of a Delaunay subgraph is the minimum 
spanning tree graph. This is tailor-made to be non-local, so that our results cannot be 
applied. 

4.1.3 Delaunay tile interactions 

In this subsection we deal with potentials on the hypergraph structure Del3 of all Delaunay 
triangles in R 2 , and we still assume that ip is bounded below and depends only on the 
hyperedge and not on the remaining configuration. Such models have been considered 
recently in [B] and [5]. We improve the existence results given there. 

4-1-3.1 Polynomially increasing triangle interactions. The first example is specified by 
assuming that the interaction of a Delaunay triangle n £ De\% is controlled by its size as 
in Proposition 14.11 

(4.3) ip(r]) < Ko + k\ 5(n) a for some constants kq > 0, k\ > and a > 0, 
where 5(rj) is the diameter of the circumcircle of n. We then have a similar result. 

Proposition 4.3. Let d = 2 and ip be a pure triangle potential on De^ which is bounded 
below and satisfies (|4.3p . Then there exists at least one Gibbs measure for ip and every 

z > {(2/V^)+2 60 )e 2K0 (ae 2 K 1 ) 1 / a /(7r 6 2 ). 

If ip is bounded above, so that k\ = 0, we recover the existence result of [SJ. 

Proof. The proof is identical to that of Proposition 14.11 One simply has to note that 
S(rj) < ((2/v^3)+2g )a when n £ Del 3 (V) for some uj £T. □ □ 

4-1.3.2 Shape- dependent triangle interactions. The shape of a Delaunay triangle can be 
captured by its angles. So let f3(n) and 7(7/) respectively denote the smallest and the 
largest interior angle of a triangle rj. We assume that ip has the form 

(4.4) pfa):=0G8fa),7fa))- 

Proposition 4.4. Let d = 2 and ip be a pure triangle potential on Dels bounded below and 
satisfying (|4.4|) . Suppose there exists some 5 > such that 

sup <K/#,7) < 00 • 

7>/3>(tt/3)-<5 

Then there exists at least one Gibbs measure for ip and every z > 0. 
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Proof. We can again apply Corollary 13.41 with the assumption (U ). The hypotheses (S) 
and (R) hold for the same reasons as in Proposition 14.11 The matrix M and the set T A 
can be defined as before, except that b is now chosen so small that /3(rj) > (vr/3) — 5 for all 
7] G Del3(a;) with w£T . Conditions (Ul ) and (U2 ) are then obvious. By Remark 13.21 
it only remains to check (XJ3 A ) because ip is scale invariant. But this is trivial. □ □ 



The conditions of Proposition 14.41 include the hard-equilaterality model, in which 
(4.5) 009,7) = 



oo if p < (vr/3) - 5, 
otherwise 

for some < 5 < 7r/3. The associated Gibbs measure produces a random Delaunay 
triangulation for which each triangle is almost equilateral in the sense that every angle 
exceeds 7r/3 — 5. This model is considered in [10J (Model 1, Section 2.3) as the crystallized 
triangulation model. 



4.2 Neighborhood-dependent hyperedge interactions 

We now turn to examples of potentials ip for which ip{r],oj) does not only depend on n 
but also on the points of to in some neighborhood of v. In two of these examples, the 
underlying hypergraph structure is based on the singleton graph 

(4.6) SG = {({x},oj) : x G u G Q} 

for which all hyperedges are singletons. Since this might seem trivial, we add immediately 
that more complex hyperedges will be hidden in the horizons of the potentials. 



4.2.1 Forced-clustering A;-nearest neighbor interactions 

In the A;-nearest neighbors model on the space M. d of any dimension d, each point x G oj 
interacts with the k points of oj that are closest to x. This model was introduced in [HE]. 
We focus here on a non-hereditary variant, in which the /c-nearest neighbors of x are forced 
to keep within distance 5 > from each other. This model was mentioned in [9] without 
precise justification of the proof, which is given here. 

Let k > 1 be some fixed positive integer. For each integer 1 < % < k, each configuration 
oj G ft with #(w) > fe+1, and each x G oj we let x t:u) be the i-th nearest neighbor of x in 
oj. More precisely, x t:ul is defined as the i-th element of oj in the total order < x on M. d , in 
which yi < x i/2 if and only if 

\yi - x\ < \y 2 - x\ or (\yt - x\ = \y 2 - x\ and y 1 < d y 2 ). 

Here, <d stands for the lexicographic order on W 1 . We also set x 0:aJ = x. 

Now let SGfc = {({x},oj) G SG : #(w) > A:+l}. We consider hyperedge potentials (p 
on SGfc of the form 



(4.7) V({x},u;) 



(x 0:w , . . . , x k ^) if diam{x 0:w , . . . , x k ^} < 5, 
00 otherwise, 



where <j) '■ (M rf ) fc+1 — > M is any bounded measurable function and 5 > a fixed constant. It is 
clear that this interaction is non-hereditary in the sense that the removal of a particle from 
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an allowed configuration can lead to a forbidden configuration. Since no kind of symmetry 
is required of (ft, the present formalism includes both the directed and undirected fc-nearest 
neighbor graphs. 

Proposition 4.5. In the forced- clustering k-nearest neighbor model described above, there 
exists at least one Gibbs measure for if and every z > 0. 

Proof. Let us apply Theorem 13.21 The hypergraph structure SG k is clearly sublinear and 
(f is bounded below. The stability assumption (S) thus follows from Remark 13.31 Next, 
for given ({x},lo) G SG^, the closed ball B(x, \x k:uJ — x\) with center x and radius \x k:uJ — x\ 
can serve as a horizon of ({x},uj). Since the corresponding open ball contains at most k 
points of u), assumption (R) follows immediately. Concerning (U), let b > be a number 
to be specified later and M = aE as in Remark 13.11 for some a > 2(6+<5). We set 

T = {uj = {x ,xi, . . . ,x k } : x G B(0,b), G B(xq,5/2) Vi = 1,. ..,&}. 

For each u G V and x G oo it then follows that \x — x k:uJ \ < 5. So (Ul) holds with rr = S, 
and (U2) follows with Cp < (/c+l^l^Hoo. Finally, we have 

e z ^U z c {T) = ^l r ({x ,..., Xk })dx ...dx k 

where i/rf is the volume of the unit ball in M. d . This shows that (U3) holds as soon as b is 
chosen large enough. □ □ 



4.2.2 Voronoi cell interactions 

Here we consider potentials that depend on the structure of the Voronoi cells. This type 
of model was introduced first by Ord; see the discussion in [21]. The spatial dimension d 
is arbitrary here. Specifically, let 

SG& = |({x},u;) G SG : \Zor 0J (x) is bounded} 

be the hypergraph structure of singletons with bounded Voronoi cells and (p be of the form 

ip({x},uj) = 0(Vor w (x)) for ({x},uj) G SG b . 

For instance, 4> might depend on the number of faces, the volume, or the surface area of 
the Voronoi cells. The necessity of passing from SG to SG^ was discussed in Example 12.21 
Each ({x},uj) G SG& has as finite horizon the bounded Voronoi flower 

(4.8) A= |J 

£eDel(w): 

The following proposition includes a result of [TJ. 

Proposition 4.6. In the Voronoi cell interaction model with bounded <f>, there exists at 
least one Gibbs measure for (p and every z > 0. 
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Proof. Let us apply Corollary 13.41 again. The stability comes from the sublinearity of 
SGfr(cj) and the lower boundedness of ip, cf. Remark 13. 31 The proof of (U' 4 ) is essentially 
the same as in Proposition 14.11 with k% = 0. Concerning the range condition (R), let 
{x} G SGf,(w) and y±, y2 any two points of the set A defined in (|4.8p . Then there exist 
6,6 G Del(w) such that y x G B(^,uj), y 2 G B(£ 2 ,oj) and x G B(€i,u) nBfe.w). The 
latter means that the union of these two balls is connected. By definition, B(£i,u>) and 
S(6, oj) contain no point of u. So, condition (R) holds with £r = 2, ur = and arbitrary 
Sr. □ □ 

In the above, the boundedness of (j> was only assumed for simplicity. In analogy to 
Propositions 14.11 and 14. 3| one can consider potentials p that are polynomially increasing 
in the diameter of the cell's flower (|4.8|) . The existence result then holds only for suffi- 
ciently large activity z. It is also straightforward to consider hard-exclusion models as in 
Propositions 14.21 and 14.41 For example, let d = 2 and 



(4.9) 0(Vor w (x)) 



if Vor w (:r) has six edges (and vertices), 
oo otherwise. 



It is then clear that Theorem 13.21 applies . and by scale invariance a Gibbsian point process 
exists for all z > 0. The typical configurations of such a Gibbs process preserve the 
topology of the triangular lattice but are less regular than those of the hard-equilaterality 
model (|4.5p for small 5. The model (|4.9p may therefore be called the randomly distorted 
triangular lattice. 



4.2.3 Adjacent Voronoi cell interactions 

Here we reconsider the potential presented in Example 12.21 which describes an interaction 
between two adjacent Voronoi cells. That is, let d = 2 and the set of hyperedges of a 
configuration oj G Q be given by 

Del2,b(w) := {{x,y} G Del2(w) : Vor w (x) and Vor aJ (y) are bounded}. 

Suppose the potential p> has the form 

ip(r],uj) = <^(Vor w (a;),Vor w (j/)) for n = {x,y} G Del 2 ,&(w). 

For instance, <p(ji,uj) can either depend on the length of the common edge or on the area 
ratio of the cells Vor^(x) and Vor UJ (y). As noticed before, the Voronoi "doubleflower" 

A= |J B(Z,u). 

£€Del(w), fnj7^0 

then serves as finite horizon of (rj,u;) G Del2,b- In this setting we have the following result 
which can be proved in the same way as Proposition 14.61 

Proposition 4.7. In the adjacent Voronoi cell interaction model in two dimensions with 
bounded <p, there exists at least one Gibbs measure for ip and every z > 0. 
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As in the case of Proposition 14.61 this example can easily be extended to the case when 
(j> is polynomially increasing in the diameter of the associated doubleflower, or when <f> 
exhibits a hard exclusion that permits the configurations in T for a suitable choice of A. 
A particular adjacent Voronoi cell interaction with hard exclusion was proposed in [10 1, 
Model 3, Section 2.3. In this model, a hard exclusion forces the cells to be neither too 
small nor too large, and a smooth contribution induces a competition between the areas 
of adjacent cells. 

4.2.4 Conclusion 

The preceding series of examples presents only a selection of possible models and could 
easily be extended. For instance, having dealt with interactions acting on single Voronoi 
cells or pairs of adjacent Voronoi cells, we could proceed to triples of Voronoi cells with a 
common point, which are indexed by Delaunay triangles, or even larger clusters of Voronoi 
cells. On the other hand, the preceding interactions can be combined (i.e., added up) to 
obtain models with a richer interaction structure. A further universe of models opens up 
if one passes to marked configurations as in Remark 13.71 Future will show which kind of 
interaction will prove suitable for modeling geometric structures in the plane or in space 
that occur in the sciences. In any case, it should be evident that the conditions of our 
theorems are flexible enough to guarantee the existence of Gibbsian point processes for a 
large variety of geometric interactions. 

5 Proofs of the theorems 

Let £, (p and z be fixed throughout this section. Before we enter into the proofs of the 
theorems, it is necessary to verify a basic ingredient of the theory of Gibbs measures, 
namely the consistency of the finite- volume Gibbs distributions. 

Lemma 5.1. Let A C A <s R d and uj G £lf ' 2 . Then 

G z A ,U n *' Z ) = 1 and J f dG ku = J (/ f dG U) G Z A,Jdu) 
for all measurable functions /:£)—>■ [0, oof. 

Proof. Let A, A and u be fixed. Since G Z A u does not depend on w/\, we can assume 
for notational convenience that ui C A c . Consider any two configurations Q G fl\ and 
£ G ^a\a- It follows straight from the definition that £a(Cu£Uw) C £a(CU£ Uw). Hence, 
sum of -£^a,£uoj(C) and a term in R U {oo} which does not depend on 

So > H A,$ju(0 < 00 when H A,lu(C U < oo, and ff A) £uw(C) < 00 when #A,w(C u < 00 • 
On the other hand, it follows that 

^a,^(C u + Ha^uUC) = h a AC u + AaWO 
for all C G Taking the negative exponential and integrating over we thus find that 
(5.1) e -HAAW = e-tfWO Zi A A0 > 
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where ^AAw(0 = / e~ ~ Ha '<"^' u ^ rR(d£') is the partial partition function for which the 
configuration £ is held fixed. Since G A w is concentrated on the set {H A w < oo, H& )U < 00} 
which is contained in {H A . Uw < 00, Sa, -uw < °°}; we can conclude that 

Ga^V = 0) = Gi iW o pr- 1 XA (^ Ai£J (-) = 0) 

= ( Z k,u)~ l j Z A,A,u(') 1 {^X,a, w (-)=0} rfII A\A = • 

Likewise, G Alu [Z a . = 00) = because J ^AAwH^^A = < 00. Combining 

these results we obtain the first assertion of the lemma. As a consequence, we can divide 
Eq. (|5.ip by Z A u Z A ^ UuJ to obtain the consistency equation G A u = J G A % G A w (dcD). □ 

We now turn to the proof of Theorem 13.21 So we assume that the hypotheses (S), 
(R) and (U) are satisfied. As usual, we construct a Gibbs measure as a limit of Gibbs 
distributions in suitable boxes. We choose M and T as in hypothesis (U) and consider for 
each n > 1 the parallelotope 

A„= (J C(k) , 

where L n = {— n, . . . , n} d . Let uJ G T be a fixed pseudo-periodic configuration with 
su Pfcez d Nc(k)(&) < 00; the last condition can be satisfied by letting U be periodic. As- 
sumption (Ul) implies that ZJ G S^-™- Combined with (|3.1|) . (S) and (U2&3), this shows 
that cJ is admissible for A n and 2; cf. (|5.8p below. So we can define the Gibbs distribution 

Gn = G z An)Z3 o pr ^ 

in A n with boundary condition ZJ and activity 0, projected to A n . Since we aim at con- 
structing a shift-invariant Gibbs measure, we will introduce a spatial averaging of G n , and 
it is convenient to work directly on the set of all shift-invariant probability measures on 

n. 

So, let P n be the probability measure on (fi, T) relative to which the configurations 
in the disjoint blocks A n + (2n+l)Mfc, k G are independent with identical distribution 
G n . We consider the averaged measure 

(5.2) P n = — [ P n o ft- 1 dx, 

where v n = \A n \ is the volume of A n . By the periodicity of P n , P n is shift-invariant. 
Moreover, the intensity i{P n ) = J N An dG n /v n of P n is finite because 

J N An e~ HA "^dn z An < e cs * {dA ^ j N An e csNA " dU z An < 00 

by (pTE]) and (S). So, P n £ &>e- 

We will show that the sequence (P n ) has an accumulation point in a suitable topology. 
As in [T3], we will take the required compactness from the compactness of the level sets 
of the specific entropy. Let us recall the necessary concepts. 

A measurable function / : O — >■ R is called local and tame if 

f(io) = f(uj A ) and l/MI < aN A (uo) + b 
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for all oj G 0, some A <g R and suitable constants a,b > 0. Let £ be the set of all 
local and tame functions. The topology of local convergence, or £-topology, on &q is then 
defined as the weak* topology induced by C, i.e., as the smallest topology for which the 
mappings P \— > J f dP (with / G C) are continuous. Note that the intensity P \— > i(P) is 
continuous in the /^-topology. 

Next, for any P G &q let P An = P ° P r A* be the projection of P to f^A n and 

up \ u z )= Uf^fdn Z A n if ^A„«n^ with density/, 
1 Anl A ' lj \ oo otherwise 

the relative entropy of P An with respect to 11^ ; here, "<C" stands for absolute continuity. 
The specific entropy of P (relative to IT) is then defined by 

(5-3) I Z (P)= lim v- 1 1(P A JW A ); 



see [12] and [15] for the existence of the limit and further properties of P . Our key tool 
is the following result of |13^ Lemma 3.4], which is based on |15t Prop. 2.6]. 

Proposition 5.2. For all c±,C2 > and z > 0, the set 

{Pe &e--I z (P)-c 1 i(P)<c 2 } 

is relatively sequentially compact in the C-topology. 

In view of this fact, the following entropy bound implies that the sequence (P n ) has a 
convergent subsequence. 

Proposition 5.3. In the limit n — )• co we have 

(5.4) P(P n ) - c s i{P n ) < \C\-\cr - lnU z c (T)) + o(l). 

Proof First of all, the definition of P n readily implies that 

(5-5) I z (P n ) = v- 1 I(G n \ni n ); 

see the proof of [12l Proposition (16.34)]. Likewise, i(P n ) = v~ x J N An dG n . By the 
definition of G n , we know further that 

(5-6) I(G n \Ul n ) = -j H hnP dG n - In . 

So we need to estimate the two terms on the right-hand side. As for the first term, 
hypotheses (S) and (Ul) give 

(5.7) J H An & dG n > -cs J N An dG n - c s , 

and the assumption on co implies that ^(d^JU) = o(v n ). 

It remains to estimate the partition function Z AntZ j. By (Ul) we can find a number 
m > 1 such that <9A^ C A n+m for all n > 1. Fix any n and let C G Q An be such that 
C ■= C U <^A= ^ r. We claim that 

(5.8) H A ^(()<c r #L n + o(v n ), 
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where the error term is uniform in £. Indeed, since n G £a„(C) when n G £(C) an d 
G fj n L n , we can write 

^mo=E E f + E E f ■ 

In view of (U2) and translation invariance, the first term on the right is not larger than 
cr #L n . Likewise, the second term is dominated by cjt ^{L n+m \L n ). This proves 
and leads us to the estimate 

(5.9) > J t T (C)e- H ^QulJdO > e - c r#^-oK) ^(r)*^. 



Combining this with (|5.5p to (|5.7p we end up with (|5.4p . □ 

The two propositions above imply that the sequence (P n ) admits a subsequence that 
converges to some P G in the £-topology. The limit P is non-degenerate, in that 
P 7^ 5$. Indeed, in view of the lower semicontinuity of I z (implied by Proposition 15. 2p 
and the continuity of the intensity i we obtain from (|5.4j) that 

(5.10) P(P) - c s i(P) < ICr^cr - lnn^(r)) . 

But hypothesis (U3) ensures that the quantity on the right-hand side is strictly less than 
z = I z (5d,) - c s i(5$). 

It is natural to expect that P is the Gibbs measure we are looking for. Unfortu- 
nately, however, we are unable to show that P is concentrated on the admissible config- 
urations. However, since P is non-degenerate, we can consider the conditioned measure 
P = P( - |{0} c ) G with P({0}) = and apply Proposition 13.11 Let us give a more 
precise statement of this proposition. 

Let £R,nji,5ji be the constants introduced in condition (R). Also, let 5- and 5 + be 
the diameters of the largest open ball in C and of the smallest closed ball containing 
C, respectively. Fix an integer m > 6£r5+/5-. For each n > 1, we decompose the 
parallelotope A n := A n+ ( 2n +i)m into the (2m+l) rf translates := A n + (2re+l)M/c of A n , 
where k G L m . For any A d W 1 let n\ > 1 be the smallest number with A„ A D A and 

> 5r/65+. For all n > n\ we consider the events 

(5.11) = { o min m iV A , > n fi } G ^ 

as well as f2„ = Un>n A ^cr' n G J-"a c - We then have the following result. 

Proposition 5.4. Given any A (s R d , we /lave O^. C O^. and dA(oS) C A n when u> G r^CT™ 
for some n > r&A- Moreover, P(f^.) = 1 /or a// P G <^e ^i/i = 0. 

Proof. Let A and n > ha be fixed and consider any uj £ Q for which the range of <p from 
A is not confined within A n . Then there exists a configuration Co € Q with cD = oj on 
A n \ A, a configuration £ G f^A and a hyperedge rj G 5a (C U w A c ) such that (p(rj, ( U wa c ) / 
ip(r),C U ^A c )- So, every horizon A of (rj,( U wa c ) as in (|2.ip hits A£. By (|2.4p . A hits 
A too. Now let A be chosen as in (R). We pick some x G A n A and y G A \ A n and, 
as in (R), a chain of I < £r balls that hit each other successively and run from x to y. 
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There is a first ball Bp hitting A£j. Shrinking Bp if necessary, we find a connected chain 
B\ , . . . , Bp of at most l R balls that are all contained in A n , connect x to dA n , and have 
either diameter at most 5 R or contain at most n R particles of £ U wa c • Since the distance 
between x and dA n is at least m(2n+l)<5_, at least one ball Bi has a diameter exceeding 
3(2n+l)5 + . As n > Sr/65 + , this bound is larger than both 5 R and 3diamA n . So there 
exist at least two indices k, k' G L m such that A^ and A 1 ^ are included in Bi and thus hold 
at most ur points of £Uwa c - At least one of these parallelotopes is different from A n , say 
A k n . This proves that u) ^ Cl cr ' n and completes the proof of the first statement. 
To prove the second claim let P G be such that -P({0}) = 0. We have 

1 - P(Cli) = p( n < n mf A E P ( N A* Z n R ) . 

n>n,A ~ O^fcGim 

By translation invariance, the last expression is equal to 

(#L m - 1) inf P(N An < n R ) . 

n>n\ 

But this term vanishes because 

P{N An < n R ) P(N Rd < n R ) = P({0}) = 

as n — > oo. The next to last identity comes from the well-known fact |19t 6.1.3] that 
P(0 < N^d < oo) = when P is translation invariant. The proof is therefore complete. □ 

The final step in the proof of Theorem 13.21 is as follows. 

Proposition 5.5. The conditional probability P = P( • |{0} c ) G is a Gibbs measure 
for £, ip and z. 

Proof. Since P G &q with P({0}) < 1, P is well-defined and belongs to To show 

that P is a Gibbs measure we fix some A <s M. d and consider the sets &cr' p defined in (|5.1ip 
for p > riA- We also set &cr~ p = Un=n A ^cr™- It is sufficient to show that 



( 5 - 12 ) L,<J dP= /~A,< P A,J AdP 

whenever / : Q, — >• [0, 1] is a local function and p > n\ is so large that / is P^ -measurable. 
Here, /a is defined by 

Indeed, letting p — > oo and setting / = 1 we then find that P(fi£.nn* '*) = P(Cl A ) and thus 



■id 



P(fi* ,z ) = 1 by Proposition For arbitrary / we obtain further that P = f G A w P(di 
Since A is arbitrary, this means that P is a Gibbs measure. 

To prove (|5.12p let / and p > n\ be fixed. It will be convenient to replace the sequence 
(P n ) introduced in (| 5 . 2 j) by an alternative sequence of measures with the same limit P. 
Suppose n is so large that A p C A n and let 

A° = {x G R d : A p + x C A„} 
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be the "A p -interior" of A n , which coincides with the closure of A n _ p _( 2 p+i)m- We define 
the (subprobability) measure 

G n '■= / G An ZJ o fi^, 1 dx = I G An _ x ^ xZ jdx ; 

v n J A° n v n Ja° 

the equality comes from the shift-invariance of ip and the symmetry of A° . The argument 
in [151 Lemma 5.7] then shows that J f dP n — J f dG n — > for all / E C. This means 
that P can also be viewed as an accumulation point of the sequence (G n ). Now let 
x E A°, so that A p C A n — x. Using the consistency lemma 15.11 and the fact that 
&cr'- P E F AAA C F( An -x)\A we find 



f 



fdG An _ x ^= f (f f dG z AuJ ) G An _ xl9xZJ (doj) , 

Jut' nQcr- p J 

and averaging over x yields 

(5.13) f fdG n = [ f A dG n . 

The integrand of the integral on the left is measurable with respect to * A and thus 
belongs to C. By (A.4) in the appendix, the integrand on the right of (|5.13p is measurable 
with respect to the universal completion F*^ ^ and thus can be squeezed between two 

functions in C which coincide P-almost surely; cf. [5], Proposition 2.2.3. So, (I5.13P gives 
(|5.12p in the limit when n runs through a subsequence for which G n tends to P in the 
/^-topology. □ 

The proof of Theorem 13.31 requires only two minor observations. First, we note that 
(til) and (U2) together with (R) imply (Ul). Indeed, let to E f. Condition (U2) then 
shows that H AmtU) (u Am ) < oo for all m > 1. If m is so large that av m — b > hr for 
the constant n/j in (R), the lower density bound (Ul) thus gives that each translate Aj^, 
k E contains more than ur points of to. Hence, every ball with at most tir points has a 
diameter no larger than 2diam(A m ). Invoking the range condition (R), we can therefore 
conclude that (Ul) holds with rr < £r max(5ft, 2 diam(A m )). 

Next we note that the non-rigidity condition in its strong form (U3) was only used 
below ()5.10p when we showed that the accumulation point P is non-degenerate; for all 
other purposes, the weak form (U3) was sufficient. However, the non-degeneracy of P is 
trivial under (Ul) because i(P n ) > a — b/v n and thus i(P) > a > by the continuity of i. 
This completes the proof of Theorem 13.31 

Appendix: Measurability 

Here we collect and prove the measurability properties we have used and add a further 
comment on measurability. Let A <s= R d be fixed. 

Claim A.l. E A := {{-q,U)) : r] E £ A (o;)} E (J/® J)*. 
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Indeed, consider the measurable functions f\(r], C> w ) = ( r ?)CU^A c ) and g{r],Q,ui) = (n,uj) 
from $7j x Oa x $7 to ttf x 17. Since £ and ip are measurable by assumption, the event 
£\ := {tp o g zfz (p o /a} then belongs to Fy ® Fa ® J 7 , and £a is equal to the projection 
image g(£\). Since .Fa is known [T71 [19] to be the Borel cr-algebra for a Polish topology 
on Qa, one can apply Prop. 8.4.4 of [5] to conclude that £\ is universally measurable, as 
claimed. 

Claim A. 2. The functions ((,uj) — > Ha,w(C) an d (C> w ) — ► H7 U (C) are measurable with 
respect to (F A (8) Fa<0*. 

To see this, we observe first that the mapping /a in Claim A.l is measurable from F/(g>F A (g> 
Fa= to F/<8>F, and therefore also from (F/(8>Fa<X>Fac)* to (F/®F)*; see |5), Lemma 8.4.6. 
In view of Claim A.l, this means that the indicator function 1^ A (?/, CU wa c ) is measurable 
with respect to (F/ (8) F A (8) Fa=)*- Given any probability measure on £lf x J1a x ^A c , 
we can therefore squeeze this indicator function between two Fj (g> F A (g) FA^-measurable 
functions which coincide almost surely; cf. Proposition 2.2.3 of Writing 

Ha,UO = X/^aC^CUwaOv^jCUwac), 

>7Cu; 

applying Theorem 5.1.2 of [19] repeatedly when the indicator function is replaced by one 
of the squeezing functions and using Proposition 2.2.3 of [5] in the converse direction we 
get the result. 

A z 

Claim A. 3. The partition function uj — > Z AuJ is measurable with respect to F Ac , £7*' € 
F A c, and G z XuJ (F) is a probability kernel from (f2^' z , F^d n A > z ) t° (fij-F). 

Let P be an arbitrary probability measure on J-\c As in Claim A. 2, we can squeeze 
the function e~ HA ^^ between two F^ <8> Fa c -measurable functions which coincide II A <S> 
P-almost surely. Integrating these functions over C with respect to II A we obtain two 
functions of uj, which squeeze Z Au} , are Fa c -measurable by the measur ability part of 
Fubini's theorem, and coincide P-almost surely. As P was arbitrary, the first result follows. 
In the same way one finds that the function uj — > n A (PT A w < oo) is F Ac -measurable. Hence 

n$>* = U z A (H Au; < oo) = 1, < Zl >w < oo} G F Ac . 

One also finds that the integral in ()2.6p depends F^c-measurably on u), which proves the 
last statement. 

Claim A. 4. Let p > n\ be fixed and suppose condition (R) holds. Claims A. 2 and A. 3 
remain valid with F A ^ A in place of Fa c as soon as all quantities are restricted to the set 

focr P defined in (|5.11|) . In particular, Q^' z D &c~ p 6 F~ , and G\ U {F) is a probability 
kernel from (n^ 2 nOcr ) ',P^ \A^ A,z nn A,p ) ^° (^)^ r )- 
Indeed, by Proposition 15.41 and (|3.1|) we have 

H A ,UO= J2 whenue^. 

The counterpart of Claim A. 2 is therefore obvious, and the analog of Claim A. 3 follows as 
before. 
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Claim A. 5. The use of universal measurability could be avoided by modifying the definition 
of £a- 



Namely, £\{oS) could be denned as the set of all rj G £[oS) for which either n n A 7^ or 

IIX<g)nx((Cl,C2) G : V^.CiUwaO ^(^UwaO) > 0. 

Then £ a G .Ff (8) T by the measurability part of Fubini's theorem, and Claims A. 2, A. 3 
and A. 4 would follow without the stars referring to universal extensions. This modified 
definition, however, is less intuitive and destroys the simple monotonicity of £\{ui) in 
A which was used in the proof of Lemma 15,11 One can still show that the required 
monotonicity holds for n^^-almost all £, but this is more involved. It is also necessary 
to redefine O^. to obtain Proposition 15.41 We therefore decided to make use of universal 
measurability. 
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